Entanglement created by spontaneously generated coherence 
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We propose a scheme able to generate on demand a steady-state entanglement between two non-degenerate 
cavity modes. The scheme relies on the interaction of the cavity modes with driven two or three-level atoms 
which act as a coupler to build entanglement between the modes. We show that in the limit of a strong driving, 
crucial for the generation of entanglement between the modes is to imbalance populations of the dressed states 
of the driven atomic transition. In the case of a three-level V-type atom, we find that a stationary entanglement 
can be created on demand by tuning the Rabi frequency of the driving field to the difference between the atomic 
transition frequencies. The resulting degeneracy of the energy levels together with the spontaneously generated 
coherence generates a steady-state entanglement between the cavity modes. It is shown that the condition for 
the maximal entanglement coincides with the collapse of the atomic system into a pure trapping state. We also 
show that the creation of entanglement depends strongly on the mutual polarization of the transition atomic 
dipole moments. 

PACS numbers: 42.50.Dv, 42.50.Gy 



I. INTRODUCTION 

The generation of continuous variable (CV) entangled light 
has attracted a significant interest due to a potential appli- 
cation in quantum information science, specifically in quan- 
tum teleportation [J], quantum telecloning |0], and quantum 
dense coding fH. Continuous variables offer the possibility 
to create entanglement deterministically and different nonlin- 
ear processes have been proposed to generate CV two-mode 
entangled beams 
down-conversion 
ing processes fill- 151 




including nondegenerate parametric 
and nondegenerate four-wave mix- 
Recently, the four-wave mixing pro- 
cess has been proposed as a potential source of narrow-band 
entangled beams, an important resource for quantum memory 
storage lfl6ll and long-distance communications Il7tl . 

Of particular interest for CV entanglement are cavity QED 
systems where entanglement between cavity modes can be 
created by coupling the modes to an atomic system or non- 
linear crystal located inside the cavity lfl8U20ll . It was shown 
that for the generation of entanglement between cavity modes, 
it is essential to create a coherence in the coupling (or entan- 
gling) system. Typical systems for entangling the modes are 
multi-level atoms or nonlinear crystals where the coherence 
can be established initially by a preparation of the atoms in 
a linear superposition of their energy states or can be created 
dynamically by a suitable driving of the atoms through four- 
wave mixing II 111 — Tl3fl or Raman-type processes Il2ll - K4ll . 

The coherence is subjected to dissipation due to the deco- 
herence process and over a long time it might be difficult to 
maintain the coherence large enough for entangling the cavity 
modes. The main source of decoherence is spontaneous emis- 
sion resulting from the interaction of the atoms with the envi- 
ronment. On a microscopic scale, the spontaneous emission 
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can be reduced or even completely eliminated, but it could be 
difficult to eliminate on a macroscopic scale where one would 
like to create entanglement using macroscopic atomic ensem- 
bles. This raises an important question of how to eliminate 
the decoherence or how to maintain a large coherence in the 
presence of the decoherence. 

In this paper, we propose a system formed by a three-level 
atom located inside a two mode cavity that can generate the 
maximal stationary entanglement between the cavity modes 
in the presence of decoherence. The atom is modelled as a V- 
type system where the dipole allowed transitions can be inde- 
pendent of each other or can be correlated through the sponta- 
neously generated coherence (SGC) J25|. The atom is driven 
by an external laser field coupled exclusively to only one of 
the atomic transitions. We use the dressed-atom approach and 
show that the effective three-level system of dressed states 
comprises a suitable medium for a non-linear coupling be- 
tween the cavity modes. We work in the strong driving limit 
which assumes that the Rabi frequency of the laser field is 
much larger than the transition damping rates and the coupling 
strengths of the cavity modes to the atomic transitions. This 
prompts us to apply the secular approximation which ignores 
the coupling of the populations of the dressed states to the co- 
herences. It is known that non-secular terms, althou gh s mall 
can have a destructive effect on coherence effects l2l[ l22ll 
or may even have constructive effects and lead to interesting 
novel features l26l - l28ll . However, we are interested in fea- 
tures created by the SGC rather than features created by the 
coherence induced by the driving field and therefore neglect 
the non-secular terms. 

We consider four scenarios, where the cavity modes cou- 
ple to the same or different atomic transitions that could be 
correlated or independent of each other. The first scenario 
represents a situation in which the atomic transitions are in- 
dependent of each other and both cavity modes couple to the 
same atomic transition that, in addition, is driven by a strong 
and in general off -resonant laser field. Physically, this system 
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behaves as a driven two-level system and the driving field oc- 
curs as a dressing field for the atoms. We demonstrate that the 
necessary and sufficient conditions for generation of the max- 
imal entanglement between the modes is to create the com- 
plete population inversion between the dressed states of the 
coupling atomic system. A population difference between the 
dressed states occurs for an off-resonant driving field. Since 
for a strong driving field there is no coherence between the 
dressed states, one could conclude that the entanglement oc- 
curs without coherence in this case. However, for a detuned 
driving field, a coherence actually occurs between the two 
bare states of the system. In other words, in the bare atom 
picture, the entanglement is created with coherence. We find 
that the maximal entanglement cannot be created in this sce- 
nario since it is not possible to create a large population differ- 
ence between the dressed states and at the same time maintain 
a strong coupling between the cavity modes mediated by the 
atom. 

In the second scenario, we include the coupling between 
the atomic transitions through the SGC, a close analog on the 
schemes of quantum-state engineering by dissipation 1291 - 13511 . 
We find that in this case, the dissipation is used to create the 
required coherence in the atomic system. The maximal sta- 
tionary entanglement can be created on demand even for the 
resonant driving field by tuning the Rabi frequency of the field 
to the difference between the atomic transition frequencies. 
As a result, the atomic system evolves into a pure trapping 
state which is an asymmetric superposition of the degenerate 
energy states. The particular pure state into which the atomic 
system evolves depends upon the ratio of the damping rates of 
the atomic transitions and the detuning of the laser frequency 
from the atomic transition frequency. The trapping effect re- 
sults in the complete population inversion between the dressed 
states of the system. In other words, the maximal steady state 
entanglement is generated when the population of the atomic 
system is trapped in a pure superposition state. 

In the third scenario, we assume that the cavity modes are 
coupled to different atomic transitions. The new feature of 
this scenario is that now the generation of entanglement is in- 
dependent of the population of the dressed states. The neces- 
sary condition for entanglement is the creation of coherence 
between the atomic transitions, the coherence that can be cre- 
ated by the SGC. 

Finally, in the fourth scenario, we consider the most general 
configuration in which each of the cavity modes is coupled 
to both atomic transitions. We show that this scenario can be 
treated as a combination of the second and third scenarios, and 
find that the generation of entanglement depends now on the 
mutual polarization of the atomic dipole moments. Depend- 
ing on whether the transition dipole moments are parallel or 
anti-parallel, the entanglement can be enhanced (reduced) by 
the constructive (destructive) interference between the atomic 
transition amplitudes. 

The paper is organized as follows. We begin in Sec. [XT] with 
a description of the proposed schemes for the generation of 
entanglement between two nondegenerate cavity modes and 
derive the master equation for the reduced density operator of 
the cavity modes. In Sec. [Ill] we study the generation and 



enhancement of entanglement between the cavity modes for 
different coupling configurations of the cavity modes to the 
atomic transitions. We are particularly interested in the role 
of the mutual polarization of the atomic dipole moments and 
the conditions for the generation of a large stationary entan- 
glement between the modes. The physical origin of entangle- 
ment between the cavity modes is explained in terms of pop- 
ulation trapping in a linear superposition of the atomic levels. 
Finally, we summarize our results in Sec. [IV] 



II. GENERAL FORMALISM 

We consider a three-level atom located inside a two-mode 
cavity. The atom is modelled as a V-type system with ground 
state 1 3), and two excited states |1) and |2) separated in fre- 
quency by Ao = UJ13 — UJ23, where W13 and W23 are atomic 
transition frequencies between states |1) -s-> |3) and |2) <-> |3), 
respectively. We shall assume that W13 > u>23 so that Ao is 
positive. This choice, of course, involves no loss of general- 
ity. The atom acts as a coupling (or entangling) medium that 
couples two non-degenerate cavity modes of frequencies uj\ 
and u-2 through the interaction of the modes with the atomic 
dipole transitions |1) -H> |3) and |2) <-> |3). In addition, the 
transition |2) <-» 1 3) is driven by a strong laser field of angu- 
lar frequency and the amplitude determined by the Rabi 
frequency 2il, as illustrated in Fig.fl] The dipole moments of 
the two allowed atomic transitions can be orthogonal or non- 
orthogonal to each other. The latter case can lead to quantum 
interference effects induced by the SGC. The cavity modes 
can simultaneously couple to one of the atomic transitions or 
to different transitions. One can also arrange a situation in 
which each of the cavity modes could couple to both of the 
atomic transitions. In this case, the coupling and the resulting 
entanglement between the modes can depend on whether the 
transition dipole moments are parallel or anti-parallel to each 
other. 




FIG. 1 : Schematic diagram of the atomic levels and one of possible 
coupling configurations of the laser and the cavity fields. A laser field 
of frequency ujl drives the |3) — > 2} transition with detuning Al 
and two non-degenerate cavity modes of frequencies oji and LO2 cou- 
ple to the driven transition with detunings Si and 82 from the laser 
frequency. 
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For an open cavity in which the atom and the cavity modes 
are coupled to the outside vacuum modes, the dynamics of the 
driven atom plus the cavity modes is conveniently described 
by the density operator p, which in a frame rotating with the 
laser frequency frequency ojl satisfies the following master 
equation (H = 1) 

^-p = ~i[H c + H a + V, p] + L c p + L a p, (1) 
at 

where 

H c = —8\a\ai + 5 2 a\a 2 (2) 

is the free Hamiltonian of the cavity modes, 

H a = (Ax, + A ) A u + A L A 22 - (A 23 + A 32 ) (3) 

is the Hamiltonian of the driven atom, 

V = (.91 01 + g 2 a 2 ) A 23 + (g 3 ai + g4,a 2 ) A 13 + H.c. (4) 

is the interaction Hamiltonian of the cavity modes with the 
atomic transitions, 

2 

L c p = K j {Zajpa^ - ajjdjp - pajoj-J (5) 
3=1 



and 



L a p = -fi[A 31 ,pAi 3 }+j 2 [A 32 ,pA 23 ] 
+ T)([A 3U pA 23 ] + [A 32 , pA 13 ]) + 



H.c. 



(6) 



are operators representing the damping of the cavity-field 
modes by cavity decay with rates ki and k 2 , and of the atomic 
transitions by spontaneous emission with rates 71 and 72. The 
parameters gi(i = 1,2,3,4) are coupling strengths of the cav- 
ity modes to the atomic transitions. We assume that in general 
the modes couple with strengths g\ and g 2 to the transition 
1 2) O 1 3) and and also can be simultaneously coupled to the 
|1) -<->• 1 3) transition with strengths g 3 and 34, respectively. 

The coefficient r/ = p^/7172 is a measure of the amount of 
coherence, the so-called SGC, induced by dissipation between 
the |1) f-> |3) and |2) <-> |3) atomic transitions. The source of 
this coherence has an obvious interpretation. Namely, sponta- 
neously emitted photon on one of the atomic transition drives 
the other transition. The degree of the coherence, measured 
by the coefficient r), depends explicitly on the mutual polariza- 
tion of the transition dipole moments with p = cos 6, where 
6 is the angle between the two dipole moments. Thus, p = 
when the transition dipole moments are orthogonal to each 
other and p attains its maximal value of p = ±1 when the 
dipole moments are parallel or anti-parallel to each other. Ob- 
viously, the SGC vanishes when p = and attains maximal 
value whenp = ±1. 

The parameter A^ = lo 23 — wl is the detuning of the laser 
frequency ljl from the atomic transition frequency oj 23 , 61 = 
bJL — u>i and 5 2 = uj 2 — ljl are detunings of the cavity modes 
bj\ and u) 2 from the laser frequency, respectively; Ay = |i) (j\ 



are the atomic transition operators between energy states \i) 
and \ j), = 1, 2, 3) of the atom. 

Since the transition |2) <H- |3) is driven by a strong, nearly 
resonant laser field, it is convenient to work in the dressed- 
state picture I36l [3~7ll . We introduce dressed states, which are 
the eigenstates of the Hamiltonian (O: 



|1> = |1>, 

|2) = sin 0| 2) -cos<£|3), 
|3) = cos0|2) +sin0|3) , 



where 



2 1 ± 1 ^ 
cos = - 
2 



Ai 

2n Q : 



(7) 



(8) 



and f2o = y A| + 40 2 is the Rabi frequency of the detuned 
field. In the dressed-state basis, the operators Ay are replaced 
by dressed-state operators Rij = \i) and the density op- 
erator of the system can be transformed to the dressed-atom 
picture by the unitary transformation 



p — cxp ( iHgt ) p exp ( —iHgt 



(9) 



where 



H a = (Ax + A )i?n + n R z - haW + ha\a 2 , (10) 

and R z — (R 22 — R 33 )/2 is the population inversion operator 
between the dressed states |2) and |3). 

Applying the unitary transformation ([9j, we find that the 
commutator part of the master equation for p contains ex- 
plicitly time dependent terms that oscillate at frequencies 5\ 
and S 2 , and the atomic dissipative part contains terms oscil- 
lating with £!o and 2£lo- in the limit of large Rabi frequency 
^0 ^ 9ii 7h me oscillating terms in the dissipative part make 
contributions of order -fi/£lo, where i = 1,2. These terms 
can be neglected in the secular approximation. The errors 
of the secular approximation are of order ji/flo and <?;/Slo- 
Thus, it is reasonable to neglect these terms on time scales 
t j^ 1 when flo 3> 5;,7i- This approximation permits 
important mathematical simplifications, and "exact" solutions 
for the steady-state density matrix elements may be obtained 
that could provide immediate insight into the physics involved 
in the problem. 

Thus, the maser equation in the dressed-atom basis and un- 
der the secular approximation simplifies to 



dt P 



V,p 



Ldp + L c p, 



(ID 



where 



V = {di[sin(20)i? z + sin 2 J R23C iOot -cos 2 0i?32C- lO °*] 
+ d 2 (sm0i? 13 e> [Ao+ ^ (Oo+Ai)lt 

-cos0i? 12 e i [ A °-5( n o- A ^] t )| +H.c. (12) 
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is the interaction of the dressed atom with the cavity modes 
with 



di 
<h 



giaie 
3 3 aie 



i8\ t 
Hit 



g 2 a 2 c 
9Aa 2 <£ 



-i&2t 

5 

-iS 2 t 



(13) 



and 



cos 2 ^[-R 2 i,/3-Ri2 
+ H.C.1 



]+H.c. 



LdP = 7i (sin 2 ^[i? 3 i , pi?i 3 ] 
+ l2 sin 2 (2<t>){[R Z: pR 
+ 72 (sinV [R 32 , pR 23 ] + cos 4 </)[i?23, pRw] + H.c.) 
+ r) sin 2 (f>({R 31 ,pR 23 } + [R 32 ,pR 13 ] +H.c.) 
+ 77ocos 2 0([i?2i,/5i?2 2 ] + [R 22 ,pRi 2 ]+B.c.) (14) 

is an operator representing the damping of the dressed-atom 
system. 

Obviously, the cavity damping term remains unchanged un- 
der the dressed-atom transformation, but the atomic dynamics 
are now determined in terms of the dressed-atom operators. 
Here, we are interested in the case of the two cavity modes 
being non-degenerated i.e., u>\ ^ u> 2 , for which the time de- 
pendence of V is quite complicated. This renders the master 
equation difficult to solve exactly, except in a special case of 
a weak coupling of the cavity modes to the atomic transitions, 
cji <C Slo- In this case, we can treat the interaction as a weak 
perturbation to the strong atom-laser interaction and find, af- 
ter tracing over the atomic variables, that the effective master 
equation for the reduced density operator of the cavity modes, 
p c = Tr aP, is of the form 



Jt Pc 



3=1 
2 

E 

3 = 1 



1^2 Aj [ojot,, 



B, 



2a j p c a <! ) - a^ajp c - p c a\aj 



^ Aj (l"',/), a , - fi, a ,a~ : - ",'/]/', 
3 = 1 
2 



■c, ■ n,"']p,"\ ■ I I.e.} 



(15) 



where Si 2 = (S 2 — Si)/ '2, Aj, Bj and Aj, Bj are the real and 
imaginary parts of complex coefficients Aj,Bj, respectively. 
The coefficients Aj and Bj have obvious interpretation as ab- 
sorption and gain rates, whereas Aj and Bj are radiative shifts 
of the cavity mode frequencies. Correspondingly, the complex 
coefficients Cj and Dj determine terms representing desired 
correlations between the cavity modes. The expressions for 
the coefficients depend strongly on the coupling configuration 
of the cavity modes to the atomic transitions and also on a 
particular choice of other parameters. The explicit analytical 
forms of the coefficients for different coupling configurations 
of the cavity modes to the atoms will be given in Sec. [TTTJ 



The master equation ( TT3T > is of a form characteristic for a 
system composed of two field modes coupled to a multi-mode 
squeezed vacuum l38ll . For this reason, to quantify entan- 
glement between the modes, we shall use the Duan's crite- 
rion ll39ll . which relates entanglement to squeezing between 
the modes. If the cavity modes were initially in a vacuum 
state, which is an example of a Gaussian state, the state of the 
modes governed by Eq. (fT3T > will remain a two-mode Gaus- 
sian state for all times t. The quantum statistics properties of 
a two-mode Gaussian state are conveniently studied in terms 
of quadrature operators of the two cavity modes 



aie 



Y, = 



V2 



(a ; V 9i -a ie - l9i ) , 1 = 1,2, (16) 



where 81 is the phase angles of the modes. If we introduce two 
operators 



-X 2 



aY 1 



1. 



-Yo 



(17) 



where a is a state-dependent real number, then, according to 
the Duan's criterion, a two-mode Gaussian state is entangled 
if and only if the sum of the variances S = (( Au) 2 ) + (( Aw) 2 ) 
satisfies the inequality 



S = 2na 2 + 2m/a 2 -4c<a 2 + — , (18) 



with a 2 = y/(2m- l)/(2n- 1), n = (a[ai) +1/2, m = 
{a\a 2 ) + 1/2, and c = |(aia2)|. Since the right-hand side of 
Eq. ( TT~8b is a positive number, we may introduce a parameter 



T = S - a 2 - 4 



(19) 



and then the condition for entanglement between the cavity 
modes is that the parameter T must be negative. 

From Eqs. ( TT8l and dT9b it is obvious that in order to calcu- 
late the parameter T, it is necessary to have available the cav- 
ity field correlation functions n, m and c. These correlation 
functions are readily found using the master equation ( fT~5b . 
from which we can derive equations of motion for the required 
correlation functions and find that they satisfy a set of coupled 
differential equations 

|(at a . ) = _ ( r j+r *) (a} aj ) 



di {aia2 



+ X]{a\<A) + Xj(aia 2 ) + 2Aj, 
= - (Ti + T 2 ) (aia 2 ) + X2(a[ai) 
+ Xi{4a 2 ) + {d + C 2 ) , (20) 



where Tj = Kj + iSi 2 — (Aj — Bj) and \j = Cj — Dj. The 
set of the differential equations d20b can be easily solved for 
arbitrary initial conditions. Since we are interested in a sta- 
tionary entanglement between the cavity modes, we analyze 
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the stability condition and find that the system is stable and 
reaches its steady-state as t — >• oo when 



Re 



r x +r 2 - V(Ti-r$)* + 4xix5 



> o. 



(21) 



The above stability condition may be simplified substantially 
for particular choices of the detunings and the Rabi frequency 
such as 61,62 3> 7i and Oq 3> 7,. 



m. ENTANGLEMENT BETWEEN CAVITY MODES 

It is clear from Eq. (TT3T > that the dynamics and entanglement 
of the cavity modes are a sensitive function of the properties of 
the driven atomic system. In order to study this dependence, 
we shall examine four scenarios of the coupling configuration 
of the cavity modes to the atomic transitions, two scenarios in 
which both modes couple to the same driven atomic transition 
and the other two in which the cavity modes are coupled to 
different transitions. A particular attention will be paid to the 
role of a specific driving of the atoms and the SGC in entan- 
gling the cavity modes. 



A. The case of both modes coupled to the driven transition 

In this section, we examine the entanglement properties of 
the cavity modes when both modes are coupled to only one of 
the atomic transitions, the laser driven transition \2) <-> |3), 
as illustrated in Fig. Q] In other words, all the fields couple to 
only one of the atomic transition. This is achieved by putting 
the coupling strengths 93 and 54 in the Hamiltonian (|4]i equal 
to zero. We shall be particularly interested in the generation 
of entanglement between the cavity modes when the coupling 
system is reduced to a simple two-level system and the role of 
the spontaneous emission in coupling of the two-level system 
to the auxiliary level Therefore, we consider separately 
two cases of orthogonal (p = 0) and non-orthogonal (p 7^ 0) 
dipole moments of the atomic transitions. When the dipole 
moments are orthogonal to each other, p = 0, and then the 
atomic transition |1) -H> |3) decouples from the driven tran- 
sition. In this case, the system reduces to that of a driven 
two-level atom. On the other hand, when the dipole moments 
are nonorthogonal, p ^ 0, and then the spontaneous emission 
on the |1) -f-s- 1 3) can influence on the two-level dynamics of 
the driven \2) <H- |3) transition. 



We start by introducing the explicit form of the coefficients 
of the master equation (15[ , which read 



Ai = g\ 



> 1 (J 1 ) S in20 + f 1 */ /^ cos 4 , 



fl{Sl)P & 22 - VoPl2 ■ 4 

Si = 9l [-7 Wi) sin 20 + /i 6 / ] / 33 2 B^ 4 ■ 

L 4 Jl2(Ol)-»7o 



/l(-£l)P22 -W21 „„„4 

7w — F~s 2 cos 1 

/12M1) -Vo 



Ci = -5152 sin 20 



F 2 {5 2 ) 



fh(^2)-rf Q 



D i = T5if2 sin 2(f> 



fl( S 2)P33 
fl2(S 2 ) ril 

/r 2 (-^)-% 2 



fl{62)P22 - V0PI2 



(22) 



where 

F^iSj) - [M 32 (^)-M 22 (^)]^ 2 -[M33(^)-M 2 3(^)]pl3 

+ [M 3 4(^)-M a4 (<y i )]^ 2 , 
F 2 (6 j ) = [M3 2 (^)-M 22 (^)]^ 2 -[M 3 3(^)-M 23 (^)]^3 

+ [M 35 (6 3 )-M2 5 (6 ] )}p s 21 , (23) 



and 



with 



f 12 (±6 j )=f 1 (±6 j )f 2 (±6 j ), j = l,2, 



(24) 



fi(±5j) = 71+72 cos^ + i ( A + - (Ai + f2 ) ± Sj , 



f 2 (±6j) = 72 ( 1 + - sin 2 20 ) + i (Oo ± 6 S ) . 



(25) 



Here, p| 2 ' P33> P12 are tne steady-state values of the atomic 
density matrix elements under the condition of ignoring the 
effect of the weak coupling between the cavity modes and 
the atom, and M mn (6j) are elements of the inverse matrix 

of u(5,y. 



/ 27! + iSj 



U(5j) 



\ 



-271 cos 2 
-271 sin 2 
Vo 

m 



cf> + iSj 



2^2 sin 4 
—272 sin 4 
Vo 
Vo 





— 272 cos' 
272 cos 4 



I 



cf> + iSj 



Vo 

—rjo cos 20 
-2r]o sin 2 
b + i5j 




'7o \ 
—770 cos 20 
— 2rjo sin 2 (f. 



b* + i5j J 



(26) 



where b = 71 + 72 sin 2 + i[A - (Q - A L )/2]. 



The remaining coefficients A2 , B 2 , C2 and D2 are obtained 
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from Eq. ( |22| | by exchanging Si with —82 and g\ with 52. 
We should point out here that in the derivation of the coeffi- 
cients d22l) . we have assumed that the states |1) and |3) are 
separated in energy by Ao + (f2rj + Ai,)/2, while the states 
|1) and 1 2} are separated in energy by Ao — (£lo — A^)/2. 
Thus, in general, the dressed states are non-degenerate. How- 
ever, by varying the Rabi frequency fin or the splitting Ao, 
one may turn the states |1) and |2) into degeneracy, whereas 
the states |1) and |3) will always remain far from resonance. 
This would happen when Ao = (fin — Al)/2. As we shall 
demonstrate in this paper, the degeneracy condition is an op- 
timal condition for entanglement between the cavity modes. 

Having defined the coefficients of the master equation for 
the case of both cavity modes coupled to the driven atomic 
transition, we now turn our attention to the possibility of gen- 
erating a stationary entanglement between the modes. In do- 
ing that we shall consider separately two cases, p = and 



1. The case of p — 

Let us first determine how much entanglement can be gen- 
erated when the atom behaves as a two-level system. The 
master equation (IT~5b can be applied to this simplified case by 
putting p = 0. Figure [2] shows the entanglement measure T 
as a function of Al for ?7o = 0, fixed detunings 81 , 82 and the 
Rabi frequency fio. The figure shows that under resonant ex- 
citation, the cavity modes are separable and become entangled 
for an off-resonant excitation. The entanglement exhibits an 
interesting behavior, in that it has two maxima which occur 
for certain nonzero values of A^, and then rapidly declines 
thereafter. A small difference 812 = —0.61 between the de- 
tunings 81 and 82 is introduced to cancel the effect of the Stark 
shifts Aj and Bj . As we see from the figure, the Stark shifts 
have a distractive effect on entanglement. 
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FIG. 2: The degree of entanglement T plotted as a function of Al 
for the case corresponding to a two-level system, 173 = (74 = and 
p = 0, with 72 = 0.02, Q = 50, Si » 8 2 = 50, m = k 2 = 
0.63, <?i = <?2 = 10 and different 812'. £12 = (solid line), 812 = 
—0.61 (dashed line). All parameters are normalized to 71. 



We would like to point out that the magnitude of the entan- 
glement is not large and there are no parameter values at which 
the entanglement could reach the optimal value T = — 1. 
Moreover, the maximal entanglement occurs at large detun- 
ings, A 7, ps ±407i, at which the driving field is weakly cou- 
pled to the atoms. We shall demonstrate in the second sce- 
nario, that the magnitude can be enhanced to its optimal value 
T = —1 by coupling the two-level system to the third level. 
To summarize, we briefly discuss the parameters characteriz- 
ing the system and the ranges of these parameters experimen- 
tally accessible. The parameters are expressed in units of the 
spontaneous emission rate 7. In the case of alkali atoms, 7 is 
of the order of 10 MHz. Driving lasers used in experiments 
are usually tunable, providing for arbitrary detuning Al, so 
that the range Al < IOO7 is easily accessible. The lasers are 
sufficiently powerful to generate Rabi frequencies up to IOO7. 



2. The case of p ^ 

We now turn to illustrate the role of the SGC on entan- 
glement creation between the cavity modes. We assume that 
the driven transition to which the cavity modes are coupled, 
is coupled by spontaneous emission to the auxiliary level 
|1). This coupling can occur for the case of non-orthogonal 
(p 7^ 0) dipole moments of the atomic transitions, and then 
the spontaneous emission on the |1) <h> |3) can influence on 
the two-level dynamics of the driven |2) <-> |3) transition. 

Since the spontaneous emission on the atomic transitions 
occurs at different frequencies and with different rates, the 
created entanglement between the cavity modes may depend 
strongly on the splitting Ao. As we shall see, the crucial for 
entanglement between the cavity modes is the relation be- 
tween fio and Ao- Figure [3] illustrates the variation of T 
with gradually increasing Ao for the case of resonant driv- 
ing, Al = 0. We see that the cavity modes become entangled 
only for p ^ and for a certain value of Ao = fin/2, the en- 
tanglement becomes optimal. In terms of the energies of the 
dressed states, the condition of Ao = fio/2 corresponds to 
the situation where the dressed states |1) becomes degenerate 
with the dressed state |2) lf40Ll4lTl . The condition of p ^ cor- 
responds to the presence of direct coupling between the states 
|1) and 1 2). Note that this coupling is induced by the dissipa- 
tive process of spontaneous emission. Since this is a resonant 
coupling, it creates a strong coherence between the states |1) 
and 1 2). Under this circumstance, the modes become strongly 
entangled and the degree of entanglement is maximal in com- 
parison with Fig. 2. The amount of the generated entanglement 
depends also on the ratio of the spontaneous emission rates, 
72/71, and the maximal entanglement of T w — 1 is achieved 
at Ao = Slo /2 and p«l for 72 -C 71 . In other words, a large 
entanglement occurs when the most of the population resides 
in the driven transition rather than in the undriven transition. 
We may summarize that by using carefully designed driving, 
such that Ao = fio/2 and carefully chosen atoms, such that 
72 <C 71, a large entanglement can be produced between the 
cavity modes via dissipation created coherence in the atoms. 

We now proceed to explain the physical origin of the pro- 
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FIG. 3: The degree of entanglement T plotted as a function of Ao 
for A L = 0, 72 = 0.02, Q = 50, 6i ~ S 2 = 50, S 12 = -0.61, 
fti = K2 = 0.63, (?i = 32 = 10, and various values of p: p — 0.98 
(solid line), p — 0.7 (dashed line), p = 0.4 (dashed-dotted line), 
p — (dotted line). All parameters are normalized to 71. 



cess responsible for entanglement of the cavity modes pre- 
dicted in the above two scenarios. As we shall see, the physics 
of the process can be quantitatively explained on the level of 
the stationary population of the atomic system. In the first in- 
stance, a simple analytical expression can be derived for the 
master equation as follows. When the frequency difference 
5 and the Rabi frequency fio are much larger than the damp- 
ing rates of the atomic transitions, S± « 82 = 8 74 and 
fio ^ 7i, the real parts of the parameters (1221) become negli- 



gible, i.e. Aj 
become Aj 



B4 



Dj m 0, and the imaginary parts 
Bj and Cj = —Dj. It is then straightforward 
to show that the master equation (TT3T > may be approximated by 



dt 



Pc = -i (Su + 2A) 



a\a\ + a\a 2 ,p c 



iD a\a\ + aia 2 , p c + L c p 



(27) 



where 



A 
D 



g 2 n (1 + cos 2 24>) 



4 (fig 



5 2 ) 



g 2 il Q sin 2(f) 

2(n 2 -5 2 ) 



[p s 2 



(P22 - P33) 

P33), 



(28) 



and, for simplicity, we have assumed equal coupling constants 
.91 = .92 = 9- 

This shows that the atomic variables contribute to the co- 
herent evolution of the cavity modes and the only relaxation 
in the system is the damping of the cavity modes. A choice of 
S12 = —2A simplifies further the master equation and leaves 
only the parametric amplifying term in its commutator part. 
This term is responsible for correlations and so for entangle- 
ment between the modes. The magnitude of entanglement at- 
tains maximal value when D maximizes. It is evident from 
Eq. (f28l > that the parameter D is different from zero only if 



the population is unequally distributed between the dressed 
states. Thus, the only one factor determines the magnitude of 
entanglement between the cavity mode, the population must 
be inverted between the dressed states of the system. For the 
case of p = 0, this can be achieved if the laser frequency is 
detuned from the atomic transition frequency 0^23- It is inter- 
esting that the entanglement is created without any coherence 
between the dressed states. There is no coherence between the 
dressed states since the Rabi frequency fio is much larger than 
all relaxation rates, fio ^> 7, , . However, we should point 
out that in the case of an off-resonat driving, there is a coher- 
ence between the bare atomic states. Thus, one can argue that 
the predicted entanglement actually occurs due to a non-zero 
coherence between the bare atomic states. 

To calculate the population inversion between the dressed 
states, we introduce density matrix elements with respect to 
the three atomic dressed states in the absence of the cavity 
modes, denoting (l|p|2) by p\%, etc. The equations of motion 
are 

Pn = — 271P11 - %(Pl2 + P21), 

P22 = 271 cos 2 cjjpn + 2j 2 (cos 4 cf>p33 - sin 4 ^22) 

+ 7] cos20(pi2 + P21), 
P33 = 271 sin 2 (j>pn - 272 (cos 4 (j>p 33 - sin 4 <pp 22 ) 



+ 2r7 sin 2 (f>(pi 2 + P21), 



P12 = — |7i + 72 sin 

- f?o(Pll + P22)- 



+ i 



A - - (fi - a l ; 



P12 
(29) 



It is evident from the above equations that the coherence P12 
induced by spontaneous emission oscillates with frequency 
Ao — (fio — Al)/2. This fact has the obvious physical mean- 
ing that the coherence attains maximal value when Ao — (fio — 
Al)/2 = 0. For Al = 0, the coherence maximizes at 
Ao = fio/2 and simultaneously the factor sin 2 4> in the co- 
efficient D equals to 1, consequently the value at which the 
entanglement, shown in Fig. [3] attains the maximal value. 

In the steady-state, the dressed state population difference 
can be worked out explicitly for both p = and p ^ 0. For 
the case of p = 0, the steady state population difference is 
given by the expression 



P22 P33 



(30) 



which clearly shows that the populations among the dressed 
states are imbalanced only for a nonzero detuning A^ 7^ 
{(f) ^ tt/4). In this case the parameter D responsible for 
the nonlinear coupling between the modes is different from 
zero. It is easy to check that the maximal entanglement seen 
in Fig.[2]is attained at the detunings corresponding to the max- 
imal value of D. Thus, we have a simple physical interpreta- 
tion of the entanglement creation by a detuned laser field. 

We stress that in the case of the detuned driving {Al 7^ 0) 
and in the limit p = 0, i.e. in the two-level situation, the 
population is unequally distributed between the dressed states, 
but it is not possible to produce atoms in a pure dressed state 
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in which |p| 2 — p| 3 1 = 1 and at the same moment having 
the coefficient D different from zero. However, for the case 
of thee-level atoms with p = 1, it is possible to have \p 22 — 
P33I = 1, in which case the population is trapped in one of 
the dressed states. The condition of the population trapping 
is unique to the SGC and can be achieved even for a resonant 
driving, A^ = 0. 

We now proceed to evaluate the population inversion when 
p = 1. A careful analysis of the steady-state solution shows 
that in the case of the level crossing at Ao = fio/2 an d m the 
limit p = 1, the population is not trapped in one of the dressed 
states but rather in one of linear superpositions 



I 1 > 




I 2 > 



I 3 > 



|,s) = a\2)+fi\l), 
\a) = /3|2) 



where 



72 sin 



7i + 72 sin 



7i 



7i + 72 sin 



(31) 



■(32) 



It is easy to check that at the level crossing condition and in 
the limit p = 1, the population is trapped in the antisymmetric 
state I a), i.e. p s aa = 1 irrespective of the detuning Al and the 
ratio between the damping rates 71 and 72. This result implies 
that the SGC is essential for the atomic system to be capable 
of achieving a pure state. In other words, the trapping effect 
is a direct manifestation of the presence of the SGC that can 
be employed to maintain the complete inversion between the 
dressed states even in the case of zero detuning between the 
laser and the atomic transition frequencies. If we incorporate 
the solution p s aa = 1 into Eq. (l28l i. we find that the resulting 
coefficient D takes the form 



D = 



3 2 ^c 



71 sin 20 



2 (fig - 5 2 



7i + 72 sin 



(33) 



from which one can easily show that the coefficient D is great- 
est when = 7r/4(Ai = 0) and 72 <C 71. This prediction 
clearly explains our numerical results presented in Fig. [3] 

To clarify the issue of the mechanism responsible for 
creation of the stationary entanglement between the cavity 
modes, we may refer to the equations of motion for the cor- 
relation functions ( l20b . It is straightforward to show that the 
limit of 6 3> Ji and f^o 3> 7i, the only damping mechanism 
of the correlation functions is the cavity damping. Thus, the 
SGC facilities correlations between the cavity modes that then 
decay with the cavity damping to a stationary entangled state. 



The case of the modes coupled to different atomic 
transitions 



We now proceed to evaluate entanglement between the cav- 
ity modes when one of the cavity modes, ai, is coupled to the 
driven |2) <-> |3) transition and the other mode a 2 is coupled 
to the undriven transition |1) -s-> |3), as illustrated in Fig. |4] 
In this case, the coupling strengths 32 = .93 = 0, then the 



FIG. 4: Schematic diagram of the coupling configuration of the cav- 
ity modes and the driven laser field. The cavity mode of frequency oji 
is coupled to the laser driven transition with detuning 61 from the 
laser frequency, while the cavity mode of frequency LJ2 is coupled to 
the undriven transition with detuning 82 from the laser frequency. 



coefficients of the master equation < fT3T > are of the form 



Ai = g\ 



P22 sm 



VoPh sm4 1 



Bi = g\ 



1 „ / „ , , sin d> 

-F 2 {61 sin 20+ ™ 9 2 
4 f2{5i)-% 



f*{-h)P22 - V0P21 cog4( 



fh(Si)-V 2 o 
C\= gigi sin <fi cos 2 tfi 

D\ = .91(74 sin cos 2 



F 3 (d2) + 

F 4 (J 2 ) - 



/r 2 (- <5 2)-?7 2 

V0P33 



/r 2 (-<52)-% 2 

with Fi (61) and i*2(5i) given in Eq. d25l l, 

F S (S 2 ) = [M 3 2(6 2 )-M 22 (S2)) p{ 2 
+ [M s5 (5 2 )-M 25 (S 2 )]p s 11 , 

F^(5 2 ) = [M 31 {S 2 )-M 21 (5 2 )]pl 2 
+ [M 35 (5 2 )-M 25 (S 2 )}p s 22 , 



(34) 



(35) 



and 

A2 = g\ 
B 2 = g\ 



h fX \ 1 M-M^k "^0^21 ■ 2 

"l (02 ) H : — ; — — n sm 



h 2 (6 2 ) 



fi2h6 2 )-v 2 
f2(-fa)ph sin 2 

fl2{-&2)-lll 



C2 = <7i<74sin0cos 



F>2 = (?i<?4sin0cos 



M<5i) + 



V0P33 



M-6i)-t$ 
h(-Si)pi 2 



(36) 
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with 



where 



hi(S 2 ) = [M 42 {-hMi + M a (-*a)Pu] cos 2 4> 
h 2 (S 2 ) = [M41 (-£ 2 )P2i + M 44 (-*a)Paa] 

fc 3 (<fl) = M43 (-<Jl)p33 ~ M 42 {Sl)p S 22 " 

ft4(<5i) = M 43 (Si)p a S3 ~ M 42 (-<Si)p| 2 " 



cos <p, 

-M45 (-<5i)Pn, 
- M 44 (-(Si)pi 2 . 

(37) 



Figure[5] shows the results for the entanglement measure T 
as a function of Ao for various values of p. Since in the case of 
p = 0, the creation of entanglement between the cavity modes 
was associated with a non-zero detuning, 7^ 0, the role of 
SGC is illustrated most clearly if one assumes a resonant laser 
field. Consequently, we choose to limit our illustration of the 
creation of entanglement to a situation in which A^ = 0. 



s_ -0.4 




FIG. 5: The degree of entanglement T as a function of Ao for the 
case of the cavity modes coupled to different atomic transitions, <?2 = 
g 3 = and 51 = g 4 = 10, with A L = 0, 72 = 2, fi = 50, <5i w 
62 = 50, <5i2 = -0.38, ki = k 2 = 0.67, and different p: p = 0.98 
(solid line), p = 0.7 (dashed line), p — 0.4 (dashed-dotted line), 
p = (dotted line). All parameters are normalized to 71 . 



As before, for the case IIII A 21 the entanglement occurs 
for p ^ and the optimal entanglement can be obtained at 
Ao = Slo/2. However, in contrast to the case A, the entan- 
glement maximizes at T w —1 for 72 = 271. It means that 
the entanglement maximizes when the transition rates of the 
dressed transition resonant with the undressed transition are 
equal. 

In order to understand this behavior of entanglement, we 
consider the coefficients of the master equation in the limit of 
5 ^> 7, and f^o ^ 7i and find that in this limit, the master 
equation (fT5b reduces to the following form 



dr 



i (tf 12 + 2A) 



a\ai + a\a 2 ,p c 



i Da\a\ + D*a 1 a 2 ,p c + L c p c , 



(38) 



A 



cos' 



Clo + 5 CL - 5 



(P22 " P33) 



D 



I 9 



sin^ (b cos^ cb 

7^_^ On - P33) + — (P22 - Pli) 

flog 2 sin 4> cos 2 A 



(Oo - 6)5 



-P12, 



(39) 



We may further simplify the master equation by choosing 
5\2 = —2 A, which leaves only the non-linear term in its com- 
mutator part. Note that comparing to the case A, there is a 
qualitative difference in the dependence of the coefficient D 
on the density matrix elements. The magnitude of D depends 
now on the coherence between the states 1 1) and |2) but not on 
the population difference. The coherence is induced by spon- 
taneous emission and can be different from zero only if p ^ 0. 
This means that the SGC is crucial for creation of entangle- 
ment between the cavity modes when the modes are coupled 
to different atomic transitions. As it is seen from Fig. [5] the 
entanglement maximizes at Ao = flo /2 and p = 1. It is easy 
to show from Eqs. d29l and OTb that for A = and 
p = 1, in the steady state the population is trapped in the anti- 
symmetric state I a). Thus, similar to the case A, the condition 
for the maximal entanglement coincides with the collapse of 
the atomic system into the pure trapping state. In this case, the 
coherence p\ 2 = —a/3 and then the parameter D reduces to 



D = - 



fl g 2 sin 2 2(j> ^7i72 
4(f7o - 6)5 7l + 72 sin 2 



(40) 



It is easily verified that the coefficient D attains its maximal 
value for (j> = ir/4 and 72 = 271. Thus, the simple formula in 
Eq. d40l > predicts accurately the parameter values of the maxi- 
mal entanglement in Fig. [5] 

In concluding this section, we would like to point out that 
the qualitative features of entanglement between the cavity 
modes depend on whether the dipole moments of the atomic 
transitions are parallel (p = 1) or anti -parallel (p = —1) to 
each other. We have already seen that in the case of parallel 
dipole moments and Ao = £Iq/2, the population is trapped in 
the antisymmetric state irrespective of the laser detuning A^ 
and the ratio between the atomic spontaneous emission rates. 
However, for the anti-parallel dipole moments, the situation is 
different. It is not difficult to show from Eqs. (|29l and OH 
that forp = — 1 and Ao = Oo/2, the steady state populations 
of the states are 



Pa 



a 



Aa 2 /3 2 , p 33 = 0, (41) 



where a and f3 are given in Eq. (32) . It is evident that in gen- 
eral the population is redistributed between the symmetric and 
antisymmetric states and only in the case of 71 = 72 sin 2 (f> 
the population is trapped in one, the symmetric superposition 
state. A consequence of this population redistribution is the 
reduction of the entanglement between the cavity modes. This 
is shown in Fig. [6] where we plot the entanglement measure 
for p = — 1 and different ratios between the atomic spon- 
taneous emission rates. For 72 7^ 271, the magnitude of 
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the entanglement is reduced and attains the maximal value 
of T = —1 for 72 = 271. This is an another demonstration 
that the maximal entanglement between the modes is achieved 
only when two correlated atomic transitions decay rates obey 
72 = 271. 




46 48 50 52 54 

FIG. 6: The degree of entanglement T plotted as a function of Ao 
for the case of anti-parallel transition dipole moments, p = — 1, with 
A L = 0, = 50, 5i « 5 2 = 50, Ki = k 2 = 0.72, and differ- 
ent 72/71: 72/71 = 0.5 (solid line), 72/71 = 1.0 (dashed line), 
72/71 = 2.0 (dashed-dotted line), 72/71 = 3.0 (dotted line). All 
parameters are normalized to 71 . 



C. Other possible couplings of the modes to the atomic 
transitions 

Finally, we briefly comment on the other possible coupling 
configurations of the cavity modes to the atomic transitions. 
The two cases discussed above predict a large entanglement at 
practically the same conditions, with only different conditions 
imposed on the damping rates of the atomic transitions. An 
another possible configuration is to couple the cavity mode u>\ 
to the undriven transition |1) <-> |3) and the mode u>2 to the 
driven transition |2) <-> |3). One can see from Fig. [4] that this 
configuration is obtained from the case B simply by replacing 
S by —5. Thus, a large entanglement could be generated in 
this configuration for the same condition as in the case B. 

The most general configuration of the coupling constants is 
the case corresponding to all of the cavity modes simultane- 
ously coupled to both atomic transitions. It is easily verified 
that this general case can be treated as a sum of two cases 
B with opposite detuning 5. By combining the two cases 
together, we find that the magnitude of the effective coeffi- 
cient D depends strongly on the sign of the parameter p. For 
p = ±1, the effective coefficient D takes the following form 

- g 2 n sin 2 20 VtT (VtT T y/TjQ .... 

2(fig-<5 2 ) 7 1 + 72 sm^0 

We see that depending on the sign of p these two coupling 
configuration can interfere constructively or destructively re- 



sulting in an enhanced or reduced effective magnitude of the 
nonlinear process. For p = 1 the configurations interfere de- 
structively such that for 71 = 72 the effective coefficient D 
vanish. On the other hand, for p = — 1 the configurations in- 
terfere constructively which results in an enhanced amplitude 
of the nonlinear process. However, the resulting magnitude of 
the effective coefficient depends strongly on the ratio 72/71 
such that D is large for 72/71 <C 1, but becomes very small, 
proportional to y7i/72 m the opposite limit of 72/71 3> 1. 
In other words, the three-level system can strongly entangle 
the cavity modes only if the spontaneous emission rate on the 
undriven transition is much larger than that of the driven tran- 
sition. 

We finish this section with a short discussion of a possibility 
to create entanglement between the cavity modes by the SGC 
in three-level atoms in the Lambda or cascade configurations. 
As we have shown, the crucial for the maximal entanglement 
is to trap the population in a pure superposition state of the 
atoms. However, it is well known that the SGC has a con- 
structive effect on trapping of the population in a pure state 
only in the V-type atoms 1251 . In the Lambda or cascade type 
atoms, the SGC has a destructive rather than constructive ef- 
fect on the trapping phenomenon li42l l43ll . 

The crucial for the entanglement is three-level atoms with 
parallel or nearly parallel dipole moments between the two 
atomic transitions. It is difficult in practice to find V-type sys- 
tems with parallel or anti-parallel dipole moments. One of 
the possibility is to use sodium dimers, which can be modeled 
as a five-level molecule in which transitions with parallel and 
anti-parallel dipole moments can be selected I44I l45ll . An al- 
ternative solution is to engineer atomic systems with parallel 
dipole moments. For example, Zhou and Swain l46ll showed 
that transitions with parallel dipole moments can be achieved 
in a three-level atom coupled to a cavity field with pre-selected 
polarization in the bad cavity limit. Agarwal ll47ll has demon- 
strated that an anisotropy in the vacuum can lead to quantum 
interference among the decay channels of close lying states. 
Another possibility is to align the dipole moments by a slow 
motion of the atoms through the medium f48j], or to apply a 
dc field to couple the upper levels of a three-level V-type atom 
with perpendicular dipole moments 14911 . 



IV. CONCLUSIONS 

We have proposed a scheme for generation on demand of 
a steady-state entanglement between two optical modes cou- 
pled to a V-type three-level atom. We have demonstrated that 
the condition for generation of the maximal entanglement be- 
tween the modes is to create the complete population inversion 
between the dressed states of the coupling atomic system. In 
the case of a two-level atom composing the entangling atomic 
system, we have shown that the sufficient condition for en- 
tanglement between the modes is to create a population dif- 
ference between dressed states of the driven atomic transition. 
However, we have found that the maximal entanglement can- 
not be created in this system since it is not possible to create 
the complete population inversion between the dressed states 
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and at the same time maintain a strong coupling between the 
cavity modes mediated by the atom. In the case of three-level 
atoms composing the entangling system, we have found that 
a stationary entanglement can be created on demand by tun- 
ing the Rabi frequency of the driving field to the difference 
between the atomic transition frequencies. The laser field 
mediates the spontaneously generated coherence between the 
atomic dipole transitions that allows to engineer the dissipa- 
tion in such a way that the atoms evolve into a pure trapping 
state. 
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